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Comment on "Phase transition in a one-dimensional Ising ferromagnet at zero 
temperature using Glauber dynamics with a synchronous updating mode" 

II Gu Yi and Beom Jun Kii^ 

BK21 Physics Research Division and Department of Physics, 
Sungkyunkwan University, Suwon 440-746, Korea 

Sznajd-Weron in [Phys. Rev. E 82, 031120 (2010)] suggested that the one-dimensional Ising 
model subject to the zero temperature synchronous Glauber dynamics exhibits a discontinuous 
phase transition. We show here instead that the phase transition is of a continuous nature and 
identify critical exponents: /3 ~ 0, !^ ~ 1, and 2 ~ 2, via a systematic finite-size scaling analysis. 
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Recently, Sznajd-Weron [1] has studied the phase tran- 
sition in a one-dimensional (ID) Ising ferromagnet at zero 
temperature using Glauber dynamics with a synchronous 
update rule. Interestingly, it has been successfully shown 
that the system exhibits a well-defined phase transition 
as the parameter Wq, the spin flipping probability for 
the null energy difference, is changed. Whereas the ab- 
sence of the phase transition at a finite temperature in the 
standard equilibrium Ising chain model is a well-known 
textbook example, the slow magnetic relaxation in the 
spin chain systems has become a hot research topic, the- 
oretically and experimentally (see references in [1|]). 

In the Glauber dynamics of the ID ferromagnetic spin 
chain at zero temperature, the energy difference AE is 
computed for the single-spin flipped configuration, and 
the spin flip is accepted at the probability 1] 



W{AE) 



1, 


for 


A£;<0, 


Wo, 


for 


A£; = 0, 


0, 


for 


A£;>o. 



(1) 



For each site of the system, the above probability is com- 
puted and each spin flip at next time step is decided. 
The synchronous update rule means that all spins in the 
system are updated at the same time, differently from 
more commonly used method of the sequential spin up- 
date rule. The use of the synchronous update rule enables 
the system to have the antiferromagnetic ordering, dif- 
ferent from the use of the sequential update rule pj . As 
an order parameter, the density p of active bonds which 
connect different spin values (up and down) is used: 



1 ^ 



(2) 



where L is the number of spins and ai{= ±1) is the 
Ising spin at the ith site in the ID chain with the pe- 
riodic boundary condition Ui+L — en applied. The den- 
sity of active bonds is especially useful in the present 
context, since it can effectively distinguish the ferromag- 
netic ordering (p = 0) and the antiferromagnetic ordering 



(p = 1). It should be noted that the system eventually 
approaches the steady state which is either fully ferro- 
magnetic pst — or fully antiferromagnetic pst — 1, and 
that the intermediate value inbetween is not possible as 
the time i — >■ oo. Due to the zero temperature nature of 
the dynamics, the ergodicity is broken and thus the time 
average and the ensemble average of pst are not equiv- 
alent to each other. The ensemble average of the order 
parameter (pst) in the present context equals the proba- 
bility of the fully antiferromagnetic stationary state, and 
thus it can change continuously across the phase transi- 
tion at Wc- We emphasize that the continuity of (pst) 
does not contradict the discreteness of pst ■ 

Figure [TJa), corresponding to Fig. 4 in 1], displays 
the ensemble-average of the steady state value of the ac- 
tive bond density (pst) as a function of the spin flipping 
probability Wq for the null energy difference. We believe 
that our simulation results and the ones presented in [l| 
are identical and both clearly indicate the existence of 
the phase transition at Wq = W^c(~ 0.5). However, in a 
sharp contrast to [l| where a discontinuous phase transi- 
tion has been concluded, in this Comment we claim that 
the system exhibits a continuous phase transition. In 
order to support our claim of a continuous phase tran- 
sition, we apply the standard technique of the finite-size 
scaling i2] of the order parameter (pst)'- 



{pst)^L-^^-'f{iWo-W,)L'^''), 



(3) 
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where f{x) is a universal scaling function, and /3 and v 
are standard critical exponents, describing the criticality 
of the order parameter and the coherence length, respec- 
tively. In Fig. [TJb), we show the scaling collapse of {pst) 
through the use of the scaling form ([3]) . It is very clearly 
shown that the phase transition at Wc ~ 0.5 is well cap- 
tured by /3 ?a 0.0 and ly w 1.0. The well-defined value 
of I' indicates that the coherence length ^ in the system 
diverges as the critical point is approached, as in a usual 
continuous phase transition [3| . 

The relaxation time r for the system to reach the 
steady state (/? = or 1) from a random initial condition 
has also been reported in Fig. 5 of [1|]. We again measure 
T in the same way as in Jl| and plot it in Fig. HJa). As 
was already observed in [ij , we also notice that (r) scales 
as L^ at Wc- This reminds us that the dynamic critical 
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FIG. 1. (Color online) (a) Ensemble average (pst) of the den- 
sity of active bonds in the steady state versus the spin flipping 
probability Wo for null energy difference, at various system 
sizes L = 80, f 60, 320, 640, and f 280. (b) Finite-size scaling 
of data points in (a) using the scaling form in Eq. Q with 
Wc = 0.5, P — 0, and i/ = 1. The existence of a continuous 
phase transition is clearly shown. Ensemble averages were 
performed over fO* random initial spin configurations. 



again the existence of the continuous phase transition at 
Wq = Wc{f^ 0.5), which is characterized by the coher- 
ence length exponent i^ « 1.0 and the dynamic critical 
exponent z « 2.0. 

In summary, different from the conclusion in |l|, we 
have clearly shown that the Ising ferromagnet chain at 
zero temperature with the synchronous update rule ex- 
hibits a continuous phase transition at Wq = Wc(~ 0.5). 
Furthermore, through the use of the standard method 
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exponent z [2] describing the divergence of the relaxation 
time is given by t '^ L^ exactly at the critical point. Ac- 
cordingly, the observation of r ~ L^ strongly indicates 
z = 2, and we can write the finite-size scaling form of the 
relaxation time as 



{T)=L^g{{Wo-W,)L'^''), 



(4) 



where g{x) is the scaling function for the relaxation time. 
The finite-size scaling form Q satisfies both (t) ~ L^ at 
Wo = Wc and (r) ~ ^^ for Wq / Wc with the prescrip- 
tion g{x) ^ x^"^ for small x since the correlation length 
follows C ~ (Wo - Wc)-". In Fig. IDJb), the scaling col- 
lapse of the relaxation time is displayed, which confirms 



FIG. 2. (Color online) (a) The ensemble averaged relaxation 
time (r) multiplied by L~^ versus the flipping probability Wo 
for various system sizes, (b) Finite-size scaling of relaxation 
time using the scaling form in Eq. ^ with Wc = 0.5, v — 1, 
and z = 2. 



of the finite-size scaling, we have identified critical expo- 
nents (3 w 0.0, ly w 1.0, and z « 2.0. 
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